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ABSTRACT 

We report wind measurements within Martian dust devils observed in plan 
view from the High Resolution Imaging Science Experiment (HiRISE) orbiting 
Mars. The central color swath of the HiRISE instrument has three separate 
charge-coupled devices (CCDs) and color filters that observe the surface in rapid 
cadence. Active features, such as dust devils, appear in motion when observed 
by this region of the instrument. Our image animations reveal clear circulatory 
motion within dust devils that is separate from their translational motion across 
the Martian surface. Both manual and automated tracking of dust devil clouds 
reveal tangential winds that approach 20-30 m s _1 in some cases. These winds 
are sufficient to induce a ~1% decrease in atmospheric pressure within the dust 
devil core relative to ambient, facilitating dust lifting by reducing the threshold 
wind speed for particle elevation. Finally, radial velocity profiles constructed 
from our automated measurements test the Rankine vortex model for dust devil 
structure. Our profiles successfully reveal the solid body rotation component in 
the interior, but fail to conclusively illuminate the profile in the outer regions of 
the vortex. One profile provides evidence for a velocity decrease as a function 
of r" 1 / 2 , instead of r -1 , suggestive of surface friction effects. However, other 
profiles do not support this observation, or do not contain enough measurements 
to produce meaningful insights. 
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Introduction 



Measurements of the winds within Martian dust devils (Balme and Greeley 2006) are 
challenging to obtain. Some knowledge of their winds stems from meteorological instru- 



mentation on the Viking landers (Hess et al. 1977) as vortices either approached or passed 



directly over the lander. Ringrose et al. (2003) analyzed the Viking data and derived vortex 



core maximum wind speeds above 40 m s" 1 for certain dust devils. However, this analysis 
carries some uncertainty because it is dependent on predictions of the dust devil's distance 
from the instrumentation, as well as assumptions regarding its radial velocity profile. Other 



studies utilizing rover-based (Metzger et al. 1999 Greeley et al. 2006, 2010) or spacecraft 



( Stanzel et al. 2008 ) observations are typically capable of reporting only the translational 



speed at which a dust devil travels across the Martian surface. Measurements of tangential 
winds within a dust devil are rare, as Cantor et al. (2006) report a single measurement of 

-1 ; 



tangential winds just below 15 m s in a dust devil observed by the Viking orbiter. 

Here we describe wind measurements within Martian dust devils using images taken 



with the HiRISE camera( McEwen et al. 2007) on board NASA's Mars Reconnaissance Or- 



biter. Although the instrument's primary objective is the observation of geologic features 
rather than dynamic meteorological phenomena, we utilize both the instrument's method 
for building color images of the surface as well as serendipitous observations of active dust 
devils to quantify their motion. 



2. Method 

The HiRISE camera utilizes a central color swath. This region features separate CCDs 
that image the surface through blue-green, red, and IR filters. During an observation, 
each CCD does not observe simultaneously; instead, a short time interval (~0.1 s) separates 
images taken through different filters. Thus, active features, such as dust devils or avalanches 



[Russell et al. 2008), observed by the color swath appear in motion when animating the 
component images (see Supplementary Material). Tracking of dust devil cloud contrast 
features between frames yields horizontal wind measurements from orbit. 

We identified 4 candidate dust devils exhibiting sufficient motion and cloud feature 
contrast. We tracked features manually by selecting a contrast feature in one image and 
identifying the same feature in a later image. We then calculated the offset distance in 
pixels between the selected points; this offset is the basis for the velocity measurement. 
Optimal manual measurements arise from dust clouds with moderate optical depth. Low 
optical depth clouds yield faint contrast features that are difficult to track, whereas high 
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optical depth clouds can yield "washed-out" features that are similarly problematic. 

Though recent advances in automated cloud feature tracking improve our measurements 
of atmospheric winds on the jovian planets (Choi et al. 2010), this method does not con- 
sistently measure dust devil motion with success. Because these automated algorithms rely 
on optimizing correlations between cropped image portions (correlation windows), the back- 
ground static Martian surface tends to dominate the comparison, frequently resulting in 
apparent velocity measurements near zero. However, we obtain some success in areas with 
more substantial dust clouds and relatively featureless background terrain. 



3. Results 



3.1. Wind vectors 



Table [T] lists the physical properties of the dust devils in our data set. Dust devils in our 
data set range from 30 to 250 m in diameter, and from 150 to 700 m in height. \Greeley et al. 



(2010) report that these diameters are typical for dust devils on the surface of Gusev crater, 



whereas Stanzel et al. (2008) describe Mars Express observations of dust devils similar in 



size to the largest dust devil in our data. Because of orbital constraints, all HiRISE images 
observe the surface at ~3 PM local Mars time. In Gusev crater, dust devil activity on 



Mars tends to diminish by late afternoon from a peak around 1 PM (Greeley et al. 2010). 



However, Mars Global Surveyor observed considerable activity occurring in its 2 PM orbit 



(Fisher et al. 2005), and Mars Express notes a similar peak in activity between 2 and 3 PM 



local time (Stanzel et al. 2008). 



Figure [T] shows two sets of wind vectors for each dust devil: manual measurements 
representing composite translational and tangential motion, and derived measurements of 
tangential motion only. All manual wind measurements shown in Figure [T] stem from all three 
possible image pairs for each dust devil. Separation of the wind vector into its translational 
and tangential components is an issue because the overall observation is short (~0.2 s), and 
tracking dust devils across the surface over timescales of minutes is not possible. Instead, 
we estimate translational velocity by calculating the overall average wind vector from the 
composite measurements. We subtract this estimate from the composite vector and exhibit 
the residual vectors in Figure [T] as the second set of wind vectors for each dust devil. Note 
that an accurate estimate using this method depends on an even spatial distribution of 
vectors from a symmetric dust devil in order for tangential motions to cancel out. Because 
some dust devils in our study better fit this criteria than others, it is important to consider 
that asymmetrical dust devils or irregularly distributed measurements will affect the velocity 



-4- 



correction in these cases. 



Overall, our composite wind measurements typically have a magnitude between 20 and 
30 m s _1 , with a maximum near 45 m s _1 . Typically, the strongest winds occur along the 
outer edge of a dust devil, regardless of its diameter. When decomposing these measure- 
ments, we note some examples of strong translational winds. First, most measurements in 
PSP_004168_1220 are typically between 20 and 30 m s -1 but likely possess a larger transla- 
tional wind component because they originate from wispy dust clouds well outside of the cen- 
tral rotating column. In addition, composite vectors for the dust devil in ESP_013199_1900 
reveal a gradation in wind magnitude from left to right and strong curl, signaling the effects 
of a considerably strong ambient wind on the inherent dust devil circulation. Our estimate 
of the translational wind speed is ~20 m s _1 . These estimates are slightly above typical 



traverse speeds (15 m s 1 and below) reported by Stanzel et al. (2008), though both Stanzel 



et al. (2008) and Greeley et al. (2010) report strong translational winds in a few instances, 
with some dust devils traveling at up to ~50 m s -1 across the Martian surface. 

In other cases, the tangential component likely dominates the measurement. The dust 
devil in PSP_009819_2130 exhibits nearly closed circulation and a relatively low average 
wind magnitude (~9 m s" 1 ). This suggests weak ambient winds, though cropping of the 
dust devil and the spatial imbalance of measurements may yield an overestimate of the 
translational wind. Finally, the dust devil captured in ESP_021925_1650 (Figure [2j displays 
substantial motion along its outer edge that produces rainbow "fringe" effects in the color 
composite image. Though imbalance in measurement location again generates overcorrected 
residual vectors, the quasi-closed circulation in the composite wind measurements (akin to 
what is seen in PSP_009819_2130) implies that ambient winds are relatively weak. This is in 



agreement with Greeley et al. (2010), who demonstrate that dust devil traverse speeds are 
> 5 m s _1 within Gusev crater. 



3.2. Radial Profiles 

Both ESP_013199_1900 and ESP_021925_1650 yielded successful measurements from our 
automated feature tracking software. The bottom row of Figure [3] shows maps displaying 
these automated results, though only a small percentage of the measurements is shown for 
clarity. Our software performs measurements at every possible pixel in an image except 
along edges, where the width of the excluded area along the edges depends on the corre- 
lation window size. The substantially larger number of measurements (typically 2 orders 
of magnitude above the number of manually tracked measurements) as a result allows con- 
struction of meaningful radial profiles. (However, the number of independent measurements 
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is less because measurements located at adjacent pixels use a base correlation window that 
is nearly indistinguishable.) 



Sinclair (1973) reported good agreement of terrestrial dust devil measurements with 



a Rankine vortex model. In a Rankine vortex, the core exhibits solid body rotation, with 
tangential velocity linearly increasing with radius. Outside of the radius with peak tan- 
gential winds, the velocity decreases as a function of r _1 . Our measurements confirm this 
region of solid body rotation, with radial profiles at several azimuths for the dust devils 
in ESP_013199_1900 and ESP_021925_1650 showing a linear increase in tangential winds. 
When defining the circulation center as the point where tangential velocities are closest to 
zero, the radius of maximum winds is ~30 ± 5 m for both dust devils. This corresponds 
to the outer visible edge of the dust devil in ESP_021925_1650, but is slightly within the 
visible edge of the dust devil in ESP_013199_1900. Unfortunately, measurements outside of 
this area of solid body rotation are unsuccessful in most cases, and we cannot completely 
illuminate the velocity structure of the dust devil. 

One such profile, however, shown in Figure [3J upper left, reveals a gradual r -1 / 2 velocity 



decrease outside of the central vortex instead of the theoretical r l . Intriguingly, Trait et al 



(2003) observed a similar profile of tangential winds within terrestrial dust devils. Both Tratt 
et al. (2003) and Balme and Greeley (2006) attribute the r -1 / 2 profile with surface friction 



causing nonconservation of angular momentum, because the densest dust clouds tend to be 
located near the surface. We stress caveats regarding this r -1 / 2 profile, as this particular 
profile is isolated to a narrow range of azimuths; other profiles do not confirm this velocity 
structure because of insufficient measurement coverage or lack of cloud tracers exterior to 



the region of solid body rotation. As laboratory studies ( Greeley et al. 2003 ) confirm the 



Rankine velocity profile, such departures from the Rankine model are likely confined to the 
surface boundary layer and are perhaps indicative of small-scale turbulent motion or vortex 
development within the central core. 



4. Error Analysis 

Typical displacements of contrast features between image frames are 5-10 pixels, with 
some above 20 pixels. The most significant source of measurement uncertainty is our imper- 
fect manual tracking technique. We estimate a 1-2 pixel measurement error accrued when 
selecting common points within contrast features. This results in a ~2-5 m s _1 uncertainty 
for each wind measurement, except for measurements from ESP_013199_1900, where lower 
spatial resolution increases the uncertainty to ~5-10 m s -1 . Measurement uncertainty from 
the automated technique, however, is below 1 pixel. 
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Because there is a slight variation in the emission angle in each color component image 
originating from the spacecraft's motion in the time interval between frames, an important 
correction to the raw measurements is compensating for parallax motion of the dust clouds, 



as Martian dust devils can potentially reach heights of ~8 km (Fisher et al. 2005). In order 
to simplify the parallax correction, we utilize non map-projected images for analysis. In 
these images, the vertical axis of the image is parallel with the spacecraft ground track, 
with the direction of spacecraft motion towards the bottom of each image. This confines the 
induced apparent motion from parallax towards the top of each image. 

We estimated dust devil heights by measuring shadow lengths and accounting for illumi- 
nation geometry. We cannot constrain the heights of the tracked dust clouds, so for parallax 
correction we assume that they are at the mid-height of the dust devil. (An exception is 
for the features tracked in PSP_009819_2130, where darker regions of the dust devil shadow 
imaged by HiRISE outside of the color swath likely correspond to the tracked dust clouds.) 
However, if the r -1 / 2 radial profile in Figure [3] signifies that most dust clouds are closer 
to the surface, the true parallax effect may be rather small, rendering our overall velocity 
measurements as a conservative lower limit. Regardless, most offset corrections are below 1 
pixel, corresponding to a velocity correction below ~2-5 m s _1 . All results discussed in this 
paper and shown in Figures [T] and [3] properly account for parallax. 

One cautionary point is that a dust devil contrast feature may appear different in 
individual frames when observed with different CCD filters. However, contrast optimization 
and image rescaling yielded no significant difference in the measurements compared with the 
default processing pipeline. Another caveat is the time-delay integration (TDI) technique 
that HiRISE employs to improve signal-to-noise. The light collection time for a given pixel 
is as much as ~11% of the time interval between frames, though this time may vary for 
different filters within the same observation. However, there is no evidence of defects in the 
automated results from smearing, and selection bias prevents affected features from inclusion 
during manual tracking. 

Finally, HiRISE is a line-scan instrument; the time interval between whole frames applies 
for lines at the same absolute ground position. Therefore, any motion of a dust devil feature 
parallel to the spacecraft ground track (and thus, perpendicular to the CCD lines) causes 
observation of said feature at a slightly shifted interval. However, the magnitude of this 
shift is two orders of magnitude below the time interval, and we ignore it for simplicity. In 
addition, HiRISE image processing applies a slight warping to images to register the color 
component images accurately and minimize spacecraft jitter. The background terrain in all 
of our image frames is well-registered, except for PSP _009819_2130, which retains residual 
jitter of ~0.2-0.5 pixels along certain areas away from the dust devil. The jitter elevates 
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the measurement uncertainty for that image, but its effect is minor relative to the dust devil 
motion after correcting for parallax. 



Discussion 



Accurate assessments of winds within dust devils are important in order to constrain 
their efficiency at dust lifting, which carries implications into the role of dust devils in the 



Martian climate and maintaining the background atmospheric haze (Newman et al. 2002 
Basu et al. 2004). Measurements indicate that the fine, entrained dust in the Martian 



atmosphere is ~2 fxm in diameter ( Tomasko et al. 1999). However, wind tunnel studies by 



Iversen and White (1982) demonstrate that boundary-layer winds necessary to move 2 fxm 



dust via shear stress exceed those observed or predicted for the Martian surface by an order 
of magnitude. 

Thus, the precise mechanism by which dust devils raise dust into the atmosphere has 
inspired several studies. Greeley et al. (2003) built a laboratory apparatus that simulated 



terrestrial and Martian dust devil vortices in order to study particle lifting threshold. They 
found that the decrease in pressure within the dust devil core facilitates the lift of particles 
(the AP effect), reducing the threshold wind speed necessary to elevate particles compared to 



boundary-layer winds acting alone. Balme and Hagermann (2006) provide theoretical argu- 



ments that the AP effect primarily operates via pressure gradients causing sudden, explosive 
release of subsurface air pockets beneath the dust, though their study could not conclusively 
determine if Martian soil would be able to degas and re-establish pressure equilibrium, pre- 
venting lifting from occurring. 

Regardless, our measurements of dust devil velocity support a significant AP, which 



could affect dust lifting into the Martian atmosphere. Greeley et al. (2003) estimate that 



lifting of 2/zm dust in atmospheric conditions similar to Mars require tangential wind speeds 
up to 20-30 m s -1 , but are likely lower on Mars due to its reduced gravity. All dust devils 
in our study exhibit areas with tangential winds above 10 m s -1 , though the prevalence of 
winds exceeding the threshold for dust lifting is variable. In addition, estimating the pressure 
perturbation of a dust devil, arguably a more important factor in the lifting of dust particles, 
is straightforward from velocity measurements. When assuming cyclostrophic balance and 
solid body rotation in the core, tangential winds of 10-20 m s -1 yield a 0.25-1% pressure 
decrease in the dust devil core relative to ambient pressure. This agrees with several studies 



of both terrestrial [Sinclair 1973), and Martian dust devils (Metzger et al. 1999 Ringrose 
et a/.||2003l). 
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Our initial search of archival HiRISE data for dust devils yielded four other dust devils 
that were unsuitable for feature tracking because of poor contrast. However, future com- 
prehensive searches for previously observed but unknown dust devils may supplement our 
knowledge of dust devil velocities. In addition, enhancements to our automated technique 
may increase the number of effective wind measurements on the current data set. Subsequent 
HiRISE observations of regions with known dust devil activity can also further expand our 
data set, though successful observations may be rare. Given the infrequency of dedicated 
meteorological instruments on surface rovers, however, an orbital imaging strategy for the 
monitoring and measurement of Martian dust devils may be an optimal plan going forward. 
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This preprint was prepared with the AAS IATeX macros v5.2. 

Image ID Latitude Longitude (east) L s Diameter (m) Height (m) 

PSP_004168_1220 -57.9° 65.5° 259.1° ~25 ~150 

PSP_009819_2130 32.8° 199.5° 120.0° -250 -650 

ESP_013199_1900 9.8° 83.1° 269.7° ~100 ~400 

ESP_021925_1650 -14.6° 175.5° 265.3° ~50 ~150 

Table 1: Physical properties of the Martian dust devils examined in this study. Note that the 
height listed for PSP_009819_2130 is for the entire dust devil, whereas we estimate the tracked 
dust clouds to be at a height < 250 m, corresponding to the densest areas of its shadow. L s is 
the Mars-Sun angle (Martian solar longitude), with L s = 0° defined as Northern Hemisphere 
vernal equinox. 
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Fig. 1. — Composite (raw) and residual (corrected for translational motion) wind measure- 
ments overlain on IRB (infrared, red, blue-green) color composite images of the dust devils 
analyzed in this work. The center of each vector corresponds to the location of the wind 
measurement. Both the color and length of the vector represent wind magnitude. 




Fig. 2.— IRB color composite of ESP_021925_1650 
fringes towards the bottom of the dust devil. 



Dust clouds in motion produce rainbow 
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Fig. 3. — Discrete radial profiles of wind measurements from ESP_013199_1900 (residual) 
and ESP_021925_1650 (non- residual) , along with automated wind measurements of both 
dust devils. Wind measurements plotted in the profile stem from the white line shown in 
the wind vector maps. The maps show only a small percentage of wind measurements in the 
data set for clarity. 



